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1 Introduction 

Among the various results in spectral analysis for selfadjoint linear partial differential 
operators, a distinguished role is played by the Weyl law, remarkable for its generality 
and simplicity — see [5J, [IE] for a detailed discussion and numerous references. Let 
us recall here a rough statement of it in the semiclassical case. Let P = p w (x, hD x ), 
< h < 1, be the semiclassical Weyl quantization on R™ of a real- valued smooth 
symbol p belonging to a suitable symbol class and satisfying an ellipticity condition at 
infinity, guaranteeing that the spectrum of P is discrete in a fixed open set Q C R. 
When [Ei,E 2 ] C Q is an interval such that vol{p~ 1 (Ej)) = 0, j = 1,2, then the 
number N(P, [Ei, E 2 ]) of eigenvalues of P in [Ei, E 2 ] satisfies 

N(P, [E U E 2 ]) = ^— (volfr- 1 ^,^])) + (1)) , h^O. (1.1) 

The leading coefficient in the Weyl law (11. ip , given by the phase space volume corre- 
sponding to the energy range [Ex,E 2 ], captures the basic physical intuition of each 
quantum state occupying a fixed volume (2nh) n in the phase space. We remark that 
the corresponding development for selfadjoint partial differential operators in the 
high energy limit has a long and distinguished tradition, starting with the works of 
H. Weyl, see [H]. 

The situation becomes quite different in the non-selfadjoint analytic case, where the 
asymptotics of the counting function for eigenvalues may no longer be governed by 
volumes of subsets of the real phase space, and the spectrum is often determined by 
the behavior of the holomorphic continuation of the symbol along suitable complex 
deformations of the real phase space. Following [3], [10], we may consider, for in- 
stance, the complex harmonic oscillator P = l/2((hD x ) 2 + ix 2 ) on L 2 (R). Here the 
range of the symbol on R 2 is the closed first quadrant, while according to the general 
results of [28], the spectrum of P is equal to {e ln ^(k + k E N}. In this case, 
it turns out that the natural phase space associated to P is given by r(e-*iR) - 
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see also [IT] for the precise spectral results in dimension one, closely related to this 
example. 

In the paper [25], it has been shown that for large and stable classes of non- self adjoint 
analytic operators in dimension two, the individual eigenvalues can be determined 
accurately in the semiclassical limit by means of a Bohr-Sommerfeld quantization 
condition, defined in terms of suitable complex Lagrangian tori close to the real 
domain. (See also [35] for the formulation of the corresponding Weyl laws.) The 
work [25] was subsequently continued in a series of papers [12] -[H], [16], [15], all of 
them concerned with the case of non-self adjoint perturbations of selfadjoint operators 
of the form, 



with the leading symbol p e (x, £) = + ieq(x,£), (x, £) G T*R 2 . Here p is real, 

so that P £= o is selfadjoint, and we assume that both p and q are analytic, with p 
elliptic outside a compact set. The spectrum of P £ near the origin is confined to a 
band of height 0(e), and to study the imaginary parts of the eigenvalues in the band, 
following the classical averaging method [10], [8], we introduce the time averages 



of q along the classical trajectories of p. 

The main focus of the works [16] and [15] was on the case when the H p -Row is 
completely integrable. The real energy surface p _1 (0) is then foliated by invariant 
Lagrangian tori, along with possibly some other more complicated invariant sets. 
Given an invariant torus A C p -1 (0), which is Diophantine (i.e. the rotation number 
of the .f/p-flow along A is poorly approximated by rational numbers), or more gen- 
erally, irrational, then the time averages (q)x along A converge to the space average 
(g)(A) of q over A, as T — > oo. When A is a torus with a rational rotation number, 
or a more general singular invariant set in the foliation of the energy surface p~ 1 (0), 
then we need to consider the whole interval Qoo(A) of limits of the flow averages 
above. 

The principal result of [16] says, somewhat roughly, that if F e R is a value such 
that F = (q)(Ao) for a single Diophantine torus A C p~ l (0), and F does not belong 
to Qoo(A) for any other invariant set A in the energy surface, then the spectrum of P £ 
can be completely determined in a rectangle [—h s /C, h s /C]+ie[F — h s /C, F + h 5 /C] 
modulo 0(h°°), where 5 is a positive exponent that can be chosen arbitrarily small, 
and e may vary in any interval of the form h K < £ < 1. The spectrum has a 
structure of a distorted lattice, with horizontal spacing ~ h and vertical spacing 
~ eh. In the work [15], we continued the analysis of the completely integrable case 
by investigating what happens when the value Fq belongs in addition to finitely many 



P E {x, hD x ) = p w (x, hD x ) + ieq w (x, hD x ), < £ < 1, 




T > 0, 



(1.2) 
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intervals Q 00(A), corresponding to rational invariant tori A. It was shown in [15] that 
the number of eigenvalues that can be created by such tori is much smaller than the 
number of eigenvalues coming from the Diophantine ones, provided that the strength 
e of the non-selfadjoint perturbation satisfies fe<e< h 2 ^ 3+s , 5 > 0, and assuming 
that F G Q OQ (A)\(q)(A), for each rational torus A. 

The purpose of the present paper is to investigate the global distribution of the 
imaginary parts of eigenvalues of P £ in an intermediate spectral band, bounded from 
above and from below by levels such as Fq, described above. In fact, when doing 
so, we shall refrain from treating the more general configurations considered in j!5j . 
with both Diophantine and rational tori present, and shall concentrate instead on 
the simpler situation of [16] , where only Diophantine invariant tori corresponding to 
F , occur. That such a study is planned by the authors was mentioned already in the 
introduction of [T6] , and here we are finally able to present the result, giving a Weyl 
type asymptotic formula for the number of eigenvalues of P e in such a spectral band. 
Roughly speaking, the main result of the present paper is as follows: let Fj G R, 
j = 1,3, F 3 < Fx, be such that Fj = (q)(Aj), where Aj C p^(0) are Diophantine 
tori, and assume that Fj does not belong to Qoo(A) when Aj 7^ A is an invariant 
set, j = 1,3. Let E 2 < < E4 be close enough to G R. Then the number of 
eigenvalues of P e in the rectangle [E 2 , E4] + ie[F 3 , Fx] is equal to 



is flow- invariant, with dQ(0) = Ax U A3. The flow-invariant sets Q(E) C p~ 1 (E), 
E 2 < E < E/x are defined similarly — see the precise statement of Theorem 12.11 
below. We may say therefore that the imaginary parts of the eigenvalues of P £ are 
distributed according to a Weyl law, expressed in terms of the long-time averages of 
the leading non-selfadjoint perturbation q along the classical flow of p. 

We would like to conclude the introduction by mentioning the work of Shnirelman [27] 
in the two-dimensional KAM-type situation, which contains the idea of exploiting 
invariant Lagrangian tori to separate the real energy surface into different invariant 
regions, in order to study the asymptotic multiplicity of the spectrum of the Lapla- 
cian. See also [5]. In our non-selfadjoint case, the idea of using invariant tori as 
barriers for the flow becomes more efficient than in the standard selfadjoint setting, 
and the main result of this work can be considered as a justification of this statement. 
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gogne in June and December of 2009. It is a great pleasure for him to thank the Insti- 
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2 Assumptions and statement of main result 

2.1 General assumptions 

We shall describe the general assumptions on our operators, which will be the same 
as in [16], [15], as well as in the earlier papers in this series. Let M stand for either 
the space R 2 or a real analytic compact manifold of dimension two. We shall let M 
denote a complexification of M, so that M = C 2 in the Euclidean case, and in the 
compact case, we let M be a Grauert tube of M — see [9] for the definition and 
further references. 
When M = R 2 , let 

p £ = P w ( x ,hD x ,e;h), < h < 1, (2.1) 

be the h-Weyl quantization on R 2 of a symbol P(x, h) (i.e. the Weyl quantiza- 
tion of P(x, e; h)), depending smoothly on e G neigh(0, [0, oo)) and taking values 
in the space of holomorphic functions of (x,£) in a tubular neighborhood of R in 
C 4 , with 

\P(x,£,e;h)\ <C>(l)m(Re (2.2) 
there. Here m > 1 is an order function on R 4 , in the sense that 

m(X) < C (X - Y) N °m(Y) J X, Y e R 4 , (2.3) 

for some fixed Co, iVo > 0. We shall assume, as we may, that m belongs to its own 
symbol class, so that m G C°°(R 4 ) and d a m = O a (m) for each a G N 4 . 

Assume furthermore that as h — > 0, 

oo 

P(x,Z,e;h)~'*rtvip j , e (x,£) (2.4) 

3=0 

in the space of holomorphic functions satisfying (12.21) in a fixed tubular neighborhood 
of R 4 . We make the basic assumption of ellipticity near infinity, 

M*,0I > ^m(Re(x,0), \(x,$\ > C, (2.5) 

for some C > 0. 
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When M is a compact manifold, we shall take P £ to be an /i-differential operator 
on M, such that for every choice of local coordinates, centered at some point of M, 
it takes the form 

Pe= J2 aa,e&h)(hD x ) a , (2.6) 

\a\ <m 

where a a ^{x; h) is a smooth function of e G neigh(0, [0, oo)) with values in the space 
of bounded holomorphic functions in a complex neighborhood of x — 0. We further 
assume that 

oo 

a a ,e( x 'i h ) ~ a ^e,j(x)h\ h-tO, (2.7) 

3=0 

in the space of such functions. The semiclassical principal symbol po tS , defined on 
T*M, takes the form 

po, £ (x,0 = ^a a , £ , (x)r, (2.8) 
if (x, £) are canonical coordinates on T*M. We make the ellipticity assumption, 

|po, £ (z,OI > ^<0 m > MeTM, \£\>C, (2.9) 

for some large C > 0. Here we assume that M has been equipped with some real ana- 
lytic Riemannian metric, so that |£| and (^) = (1 + l^] 2 ) 1 ^ 2 are well-defined. We shall 
consider the operator P e as an unbounded operator: L 2 (M , (i(dx)) —> L 2 (M , (J>(dx)) , 
with the domain H m (M), the standard Sobolev space of order m. Here n(dx) is the 
Riemannian volume element on M. 

Back in the Euclidean case, the ellipticity assumption (12. 5p implies that, for h > 
small enough and when equipped with the domain 

H(m) := (m w (x, hD))' 1 (L 2 (R 2 )) , (2.10) 

the operator P e becomes closed and densely defined on L 2 (R 2 ). We shall furthermore 
make the following assumption concerning the growth of the order function m(X) 
at infinity: for some k G N, we have 

— k G ^(R 4 )- (2-11) 

mr 

Since P e can be replaced by its k-th power P £ fc in what follows, we may and will 
assume henceforth that k — 1. The assumption that 

— G L X (R 4 ) (2.12) 
m 

will prove very useful in various trace class considerations throughout the paper. In 
the compact case, it amounts to assuming that the order m of P £ satisfies m > 2. 
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We shall write p £ for p e and simply p for p ,o- We make the assumption that 

P £=0 is formally selfadjoint. (2-13) 



It follows from the assumptions (12.5p . (12. 9p . and (I2.13P that the spectrum of P e 
in a fixed neighborhood of G C is discrete, when < h < ho, < e < sq, with 
ho > 0, eo > sufficiently small. Moreover, if z G neigh(0, C) is an eigenvalue of P £ 
then Im2 = 0{e). 

We furthermore assume that the real energy surface p _1 (0) D T*M is connected 
and that 

dp^O along p _1 (0) PI T*M. (2.14) 

We then have the following general result, stating that the distribution of the real 
parts the eigenvalues of P e near G C is governed by the usual Weyl law of the form 
(11.11) . with the following remainder estimate, 

#{z G Spec(P e ); E 1 < Rez < E 2 } = — ±— [ [ dxd^ + O (max(e,h)h- n ) . 

(2nh) n JJ p -i i[EuE2]) 

(2.15) 

Here Ej G neigh(0,R), j = 1,2, are close enough to but independent of h. In 
fact, in (I2.15P we may allow Ej to depend on h with E 2 — E\ > cmax(e, h), for some 
c > fixed. The Weyl law (I2.15P is obtained by following the general arguments of 
section 5 of [33], which are based on an adaptation to the semiclassical case of the 
ideas in [22], [23]. See also [1]. 

Let H p = p'^-d x —p' x -d^ be the Hamilton field of p. Throughout this paper, we shall 
work under the assumption that the if p -flow is completely integrable. Following [15] , 
let us now proceed to discuss the precise assumptions on the geometry of the energy 
surface £> _1 (0) D T*M in this case. 

2.2 Assumptions related to the complete integrability 

We assume that there exists an analytic real valued function / defined near p~ l (0) D 
T*M such that H p f = 0, with the differentials df and dp being linearly independent 
on an open and dense set C neigh(p _1 (0) D T*M,T*M). For each E G neigh(0,R), 
the level sets 

A^ = /" 1 (/x)n P - 1 (E)nrM 

are invariant under the if p -flow and form a singular foliation of the three-dimensional 
hypersurface p~ l (E) (~)T*M. When (fi, E) G R 2 is such that df A dp 7^ along A^e, 
then A^e is a two-dimensional real analytic Lagrangian sub manifold of T*M, which 
is a finite union of tori. In what follows we shall use the word "leaf and notation A for 
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a connected component of some A^ E . Let J be the set of all leaves in p 1 (0) D T*M. 
We have a disjoint union decomposition 

p -1 (0) HT*M = (J A, (2.16) 

AeJ 

where A are compact flow-invariant connected sets. 

The set J has a natural structure of a graph whose edges correspond to fami- 
lies of regular leaves and the set S of vertices is composed of singular leaves. The 
union of edges J\S possesses a natural real analytic structure and the corresponding 
Lagrangian tori depend analytically on A £ J\S with respect to that structure. 

For simplicity, we shall assume that / is a Morse-Bott function restricted to 
p _1 (0) H T*M. In this case, the structure of the singular leaves is known [2], [39] , 
and the set J is a finite connected graph. We shall identify each edge of J analytically 
with a real bounded interval and this determines a distance on J in the natural way. 
The following continuity property holds, 

For every A £ J and every s > 0, there exists 5 > 0, such that if (2.17) 
A £ J, dist j(A, A ) < 5, then A C {p £ n T*M; dist(p, A ) < e}. 

By the Arnold-Mineur-Liouville theorem [I], each torus A £ J\S carries real 
analytic coordinates x\, x 2 identifying A with T 2 := R 2 /27rZ 2 , so that along A, we 
have 

H p = ai9 xi + a 2 d X2 , (2.18) 
where a\, a 2 £ R. The rotation number is defined as the ratio 

w(A) = [ai : a 2 ] £ RP 1 , (2.19) 

and it depends analytically on A £ jyS 1 . 

In what follows we shall write 

Vs =p + ieq + 0(e 2 ), (2.20) 

in a neighborhood of p~ 1 (0) PI T*M, and for simplicity we shall assume throughout 
this paper that q is real valued on the real domain. (In the general case, we should 
simply replace q below by Re q.) For each torus A £ J\S, we define the torus average 
(g)(A) obtained by integrating q\\ with respect to the natural smooth measure on 
A, and assume that the analytic function J\S 3 A i— > (<?)(A) is not identically 
constant on any open edge. Here the integration measure used in the definition of 
(g)(A) comes from the diffeomorphism between A and T 2 , given by the action- angle 
coordinates. 

We introduce 

(q) T =^J qoexp{tH p )dt, T > 0, (2.21) 
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and consider the compact intervals Q 00(A) C R, A G J, defined as in [16], [T5] . 



Qoo(A) 



lim inf (q)r, hm sup(gW 

T->oo A T->oo A 



(2.22) 



Notice that when A G J\S and a; (A) ^ Q then Q 00(A) = {(g) (A)}. In the rational 
case, we write w(A) = — , where m G Z and n G N are relatively prime, and where 
we may assume that m = 0(n). When k(u(A)) := \m\ + \n\ is the height of 00(A), 
we recall from Proposition 7.1 in [16] that 



Remark. As jyS 1 9 A — > A G S, the set of all accumulation points of (5) (A) is 
contained in the interval Qoo(Ao). Indeed, when A G J\<5 and T > 0, there exists 
P = Pt,a G A such that (9) (A) = (g)r(p)- Therefore, each accumulation point of 
(g)(A) as A — )• A G S 1 , belongs to [infA (g)T) sup A (g)^]. The conclusion follows if 
we let T — > 00. 

2.3 Statement of the main result 

From Theorem 7.6 in [16] we recall that 



as s, h, 5 — 0. Let us also recall from [16] that a torus A G J\S is said to be 
Diophantine if representing H p \\ = a±d Xl + a 2 d X2 , as in (12.181) . we have 




(2.23) 



-Im (Spec(P e )n{z;|Rez| < 6}) C inf |J Q 0O (A) - o(l), sup [j Q 00 (A) + o(l) 




(2.24) 



a • fc| > 



1 



(2.25) 



for some fixed C , A^q > 0. 
Let 



F,G|Jg 0O (A), j = 1,3, F 3 <F 1 , 



AeJ 



be such that there exist finitely many Lagrangian tori 



A 1;j ,... ,A Lj;j G J\S, J = l,3, 



(2.26) 



that are uniformly Diophantine as in (I2.25p . and such that 



(q)(A kJ )=Fj foTl<k<Lj, j = 1, 3, 



(2.27) 
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with 



d A {q){A ktj )^0, l<k<L 31 j = 1,3. 



(2.28) 



We shall make the following global assumptions, for j = 1,3: 

F 3 i |J Qoo(A). 

AeJ\{Aij,...,A t j} 



(2.29) 



Here we recall from [TB] that the earlier assumptions imply that for j = 1,3, Fj ^ 
g oo (A) for A kJ j£ A e neigh (A fcJ , J), 1 < k < Lj. 

For notational simplicity only, throughout the following discussion we shall assume 
that L\ = L3 = 1 and we shall then write Ai := A^i and A 3 := A 13 for the 
corresponding Diophantine tori C p~ 1 (0) D T*M. 

Assume that the strength of the non-selfadjoint perturbation e satisfies 

h K <e<h s , 5>0, (2.30) 



for some fixed K > 1. In the work [16], under the assumptions (I2.29p . (I2.30p . for 
h > small enough, we have determined all eigenvalues of P £ mod 0(h°°), in a 
spectral window of the form 



c" c 



+ is 



e s „ e 5 
p. j 

C" 3 C 



J = 1,3. 



(2.31) 



Here 5 > is sufficiently small and the constant C > is sufficiently large. Recall 
also that the eigenvalues of P e in the region (I2.3ip form a distorted lattice, and their 
total number is 

h 2 ' 

Our purpose here is to study the distribution of the imaginary parts of the eigenvalues 
of P e in suitable "large" sub-bands of the entire spectral band 

{z e neigh(0,C); lmz = 0(e)}. 

Specifically, we shall be interested in counting the number of eigenvalues of P £ in a 
region of the form 



£ 1 f! 

C" c 



+ ie [F 3 , Fi] 



(2.32) 



bounded from above and from below by the Diophantine levels Fj, j = 1,3, intro- 
duced above. 



The following is the main result of this work. 
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Theorem 2.1 Let P s satisfy the assumptions made in subsections 2.1 and 2.2, and 
in particular, ( I2.12p . Let Fj G UAejQoo(A), j = 1,3, F 3 < F 1; be such that the 
assumptions ( I2.27p . (I2.28P . (I2.29P are satisfied. Let 0<S<l<K<oo and assume 
that h K < e < h 5 . Let E 2 < < F 4 satisfy \Ej\ ~ e 5 , where < 5 < 1/3 is small 
enough, so that 

h < e loI \og-. 

e 

Then the number of eigenvalues of P £ in the rectangle 

R=(E 2 ,E i ) + ie(F 3 ,F 1 ), (2.33) 
counted with their algebraic multiplicities, is equal to 

^E 2 <E<E 4 



vol f |J Q(E) j + O ^£ 35 log ^ J . (2.34) 



(2nhy 
Here the set 

tt(0) = { P G^ 1 (0)nT*M ; g oo (A(p)) C [F 3 ,F!]}, (2.35) 

is flow-invariant, whose boundary is AxUA^. We define the set Q(E) C p~ l (E)(lT*M 
to be the set which is close to Q(0) with boundary Ai(F) U A 3 (F) where Aj(E) is the 
unique invariant torus C p~ l (E) DT*M close to Aj, determined, thanks to (12.281) . 
by the condition (q)(Aj(E)) = Fj, j — 1,3. 

Remark. Combining Theorem 12. II together with the Weyl law f)2.15p . we see that as 
h — > 0, the ratio 

#(Spec(P £ )ni?) 
#(Spec(P £ )n([F 2 ,£ 4 ] + *R)) 
converges to the expression 

/ p -i(o)ln(o)^feO) 
J P -i(0) £ ( rf ( x 'O) 

Here £(d(x,£)) stands for the Liouville measure on p _1 (0). 

The result of Theorem 12.11 admits a natural extension to the case when the real 
energy E e [E 2 , F 4 ] varies in a sufficiently small but /i-independent neighborhood of 
G R, which we now proceed to describe. When doing so, we shall assume that the 
tori Aj, j = 1,3, introduced in (12.261) . satisfy the following additional isoenergetic 
assumption, 

d A=A .u(A) ^ 0, j = l,3. (2.36) 
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Here a; (A) is the rotation number of A, introduced in ( I2.19p . Let us represent Aj ~ 
{£ = 0} C T*T 2 X using the action-angle coordinates near Aj, so that p = p(£), 
u (0 = '■ %>p(0]- It follows from ( I2.36P that the map 

neigh(0,R 2 ) 3 £ (p(0,<40) e R x RP 1 

is a local diffeomorphism. There exists therefore an analytic family of Lagrangian 
ton Aj(E) C p^iE) n T*M, E G neigh(0,R), j = 1,3, close to A,-, such that 
u)(Aj(E)) = u(Aj), j = 1, 3. Let us set when E G neigh(0, R), 

F j (E) = (q)(A J (E)), j = 1,3, (2.37) 

and notice that an application of Theorem 12. II together with the results of P3] allow 
us to conclude that the number of eigenvalues of P £ in the region 

Rez G [E 2 ,E 4 ]; F 3 (Rez) < <Fx(Rez)>, E 2 < < E 4 , \Ej\ ~ e s , 

(2.38) 

is given by 



( U 



vol II +° e ■ ( 2 - 39 ) 



Here fi(0) = Q(0) in fl235|) and the set C p~ l (E) n T*M is close to fi(0), with 

the boundary Ai(E) UA 3 (E). Covering a sufficiently small but fixed open interval 
J C R containing G R by (9(e~ <5 ) subintervals of length 0(e s ) and applying the 
result of Theorem 12.11 in the form (I2.39P uniformly as E G J varies, we get the 
following conclusion, by summing the individual contributions from the subintervals. 

Theorem 2.2 Let us keep all the assumptions of Theorem 12.11 and assume in ad- 
dition that the isoenergetic condition (I2.36j) holds. Let C > be sufficiently large 
and let E 2 < < E± be independent of h with \Ej\ < 1/C, j = 2,4. Introduce the 
functions 

F,(E) = (q)(A j (E)) 1 E G neigh(0,R), 

for j = 1, 3 ; where Aj(E) C p~ l (E) C\T*M are Lagrangian tori close to Aj, such that 
u)(Aj(E)) = u(Aj), j = 1,3. Then the number of eigenvalues of P £ in the region 

Im z 

E 2 <Rez< E 4 , F 3 (Rez) < —— < F x (Re z) 
counted with their algebraic multiplicities, is equal to 

(^( voi Ly< E ," (£, ) +o(£i) )- (2 - 4o) 

Here < 5 < 1 satisfies the same smallness condition as in Theorem \l.\\ 
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In the course of the proof of Theorem 12. 1[ we shall assume, as we may, that the resol- 
vent (z — Pe)" 1 exists when z G 7 := OR. We equip 7 with the positive orientation, 
and decompose 



where 7, for j = 1,3, is the part of 7, where Iraz = sFj. Correspondingly, we have 
Re z = Ej along 7^, when j = 2, 4. We are going to be concerned with the asymptotic 
behavior of the trace of the spectral projection of P £ associated to the spectrum of 



giving the number of eigenvalues of P £ in the rectangle R. 

The plan of the paper is as follows. In Section [3], we recall some basic facts and 
estimates, established in [16], to be used in the proof of the main result. In par- 
ticular, we recall the main features of the averaging procedure along the H p — flow, 
effectively replacing q in (I2.20p by (q)r in (12.211) . as well as the associated suit- 
able IR- manifold C T*M, playing the role of the real phase space throughout the 
proofs. Section H] is devoted to the analysis of the trace integrals along the verti- 
cal segments 72, 74, contributing to the trace in (12.421) . Following the ideas and 
methods of [22], [23], [32], [33], [31], here we make use of auxiliary trace class per- 
turbations, constructed so that the perturbed non-selfadjoint operator has gaps in 
the spectrum. Let us also notice that the analysis of Section H] is of a general nature 
and does not depend directly on the paper [16]. It is in Section [51 concerned with 
the resolvent integrals along the horizontal segments 71, 73, that the results of [16] 
become important. The trace analysis along the horizontal segments proceeds by 
means of a pseudodifferential partition of unity, and in particular, when understand- 
ing the contributions coming from small neighborhoods of the Diophantine tori Aj, 
j = 1,3, we apply the quantum Birkhoff normal form for P £ , constructed in [16J. 
Section [5] is concluded by combining the contributions from the different boundary 
segments to derive a leading term for the trace integral (I2.42p . The Weyl law (12.341) 
in Theorem 12.11 then follows when we let the averaging time T become sufficiently 
large. The brief Section [6] is concerned with the task of deriving an analog of The- 
orem 12.11 in the case when the Hamilton flow of p is periodic on energy surfaces 
p~ l (E) fl T*M, E G neigh(0,R), which was the main dynamical assumption in the 
series of works [12] [E]. A classical Hamiltonian with a periodic flow can be con- 
sidered as a degenerate case of a completely integrable symbol, and thus, it seemed 
natural to include this discussion here. In the final Section we apply Theorems 
12.11 and 12.21 to a complex perturbation of the semiclassical Laplacian on an analytic 
simple surface of revolution, and complement the corresponding discussion in [16] . 



7 = 71 U 72 U 73 U 74, 



(2.41) 



P £ m R, 




(2.42) 
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3 Review of some results from 
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As alluded to in Section [2J the analysis of the trace of the spectral projection in 
(I2.42p will proceed by working in the exponentially weighted /i-dependent Hilbert 
space H(A), constructed and exploited in [16]. The purpose of this section is there- 
fore to recall briefly the definition and the main features of the weighted space 
H(A), following [IB], as well as some resolvent estimates for P e , when viewed as an 
unbounded operator on H(A). 

Following the discussion in Section 5 of [16], let us recall the following result. 

Proposition 3.1 Let us keep all the general assumptions of Section® and in par- 
ticular (T2J27D, (|2T28]1 . and fl2^9l) . Let 

Kj : neigh(A j; T*M) neigh(£ = 0, T*T 2 ), j = 1, 3, (3.43) 

be the real and analytic canonical transformation given by the action-angle coor- 
dinates near Aj, such that Kj(Aj) = {£ = 0}. Then the leading symbol p e of P £ 
expressed in terms of the coordinates x and £ on the torus side, takes the form 

Pj (0+ie qj (x,0 + O(e 2 ), Pj (0 = a r £ + O(e), J = 1,3. 

Define 

(<?j)(0 = j^y J qj{x,0 dx, £ G neigh(0,R 2 ), 

so that (qj)(0) = Fj, and d P j(0) = a,j and d(qj)(0) are linearly independent, for 
j = 1,3. Let < e <C 1 be such that e ^> max(e,h). Then there exists a globally 
defined smooth IR-manifold A C T*M and smooth Lagrangian tori A\ and A3 C A, 
such that when p G A is away from an ^-neighborhood of Ai U A 3 in A then 

|ReP E (p)|>^- (3.44) 
or _ _ 

EE EE 

llmP^-eF,]^—- and |ImP e (p) - eF 3 \ > —. (3.45) 

The manifold A is O(e) -close to T*M and agrees with it away from a neighborhood 
0/^(0) flT*M. We have 

P £ = 0(1) : H(A,m) -> H(A). 

When j = 1, 3, there exists an elliptic h-Fourier integral operator 

U 3 = 0(1) : H(A) L 2 (T 2 ), 



14 



such that microlocally near Aj, we have 

U,P £ = (pf>{hD x , e; h) + R N+1J (x, hD x , e; fr)) U r 

Here pj (hD x ,e;h) + Rj^ + ij(x,hD x ,s;h) is defined microlocally near £ = in 
T*T 2 ; the full symbol of pj (hD x , e; h) is independent of x, and 

R N+1)j (x,£,e;h) = 0((h,£,e) N+1 ). 

Here the integer N is arbitrarily large but fixed. The leading symbol of Pj (hD x , e; h) 
is of the form 

Remark. In the work [16], the case of a single Diophantine level Fj has been con- 
sidered, with the definition of the Hilbert space H(A) depending on Fj. It is clear 
however, from the arguments of [IE] , that the construction of the IR- manifold A can 
be carried out so that it works for both F\ and F$, since the corresponding Diophan- 
tine tori Aj, j = 1, 3, are disjoint. Let us also observe that the tori Aj C A are 
C(e)-close to A j C p -1 (0) n T*M in the C°°-sense, j = 1, 3. 

Remark. Let us recall from [TH] that the space L^(T 2 ) of Floquet periodic functions, 
occurring in the statement of Proposition 2.1, consists of all u e Lj 2 oc (R 2 ) such that 

Q 1 

u(x -v) = e ie - u u(x), 6 = —- + ^, v G 2ttZ 2 . 

Here S = (Si,^) is the pair of the classical actions computed along two suitable 
fundamental cycles in Aj, and the 2-tuple ko E Z stands for the Maslov indices of 
the corresponding cycles. 

Let Gt be an analytic function defined in a neig hborhood of p _1 (0) n T*M, such 
that 

Hp G T = q - (q) T , (3.46) 
where T > is large enough but fixed. As in [16J, we solve (I3.46P by setting 

G T = j TJ T (-t)q o exp (tH p ) dt, J T (t) = f^j , 

where the function J is compactly supported, piecewise linear, with 

J'(t)=5(t)-l hm (t). 
It follows from the results of [16] that there exists a C°° canonical transformation 
k : neigh(p _1 (0), T*M) -> neigh(p _1 (0), A), (3.47) 
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such that, in the case when M = R 2 , 



k( P ) = p + ieH G (p) + 0(e 2 ), (3.48) 

in the C°°-sense. Here the function G G C^°(T*M) is such that in a neighborhood 
of p _1 (0) fl T*M, away from a small but fixed neighborhood of Ai U A3, we have 
G = Gt- Near Aj, j = 1,3, the function G is of the form G o kJ 1 , where Kj is 
the canonical transformation near Aj given by the action-angle variables, defined in 
(I3.43p . and G is an analytic function defined in a fixed neighborhood of £ = 0, such 
that 

H p G = q-r, r(x,0 = (g)(0 + ^ 7V ), 

where JV G N is arbitrarily large but fixed. We refer to Section 2 of [TBj for the 
details of the construction of the weight function G, also in the compact case. 

Remark. From [16], we know that in a complex neighborhood of p _1 (0) fl T*M, away 
from a small neighborhood of Ai U A3, we have 

A = exp (ieH GT ) (T*M) . 



We now come to recall the definition of the Hilbert space H(A). When doing so, 
let us first concentrate on the case when M = R 2 . We shall then take the standard 
FBI-Bargmann transform 

Tu(x) = Ch~ 3/2 J e~^ x - y)2 u{y) dy, C > 0, (3.49) 

and remark that according to [50] , for a suitable choice of C > in f l3.49p . T maps 
L 2 (R 2 ) unitarily onto 

if^C 2 ) := Hol(C 2 ) n L 2 (C 2 ; e~^L{dx)). 

Here $o(x) = \ (Imi) 2 and L(dx) is the Lebesgue measure in C 2 . 

When viewing T in ( I3.49P as a Fourier integral operator with a complex quadratic 
phase, we introduce the associated complex linear canonical transformation 

kt ■ (y,v) ^ (z,0 = (y - 177,77), (3.50) 
mapping the real phase space R 4 onto the linear IR-manifold 



MA) = {|i,-? £ ];ieC I |=:A,, (3.51) 



It was shown in [16] that the representation 

i dx 
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holds. Here the function $ e G C°°(C 2 ;R) is uniformly strictly plurisubharmonic, 
with <3> £ — $ compactly supported, and such that uniformly on C 2 , 

$ e (x) - $ (z) = 0(e), V ($ e (x) - $ (x)) = O(e). 

The Hilbert space H(A) is then defined so that it agrees with L 2 (R 2 ) as a linear 
space, and it is equipped with the norm 

IMI = IMkA):=H^IU| e , L% a = L 2 (c 2 ;e-^L(dx)^ . (3.52) 

Furthermore, the natural Sobolev space H(A,m) C H(A), associated to A and the 
order function m is introduced so H(A,m) agrees with the space H(m) in ( 12.101) as 
a space, and it is equipped with the norm 

\\ u \\ 2 H(A,m) = J \Tu(x)\ 2 m 2 (x)e- 2 -^ £l L(dx), (3.53) 

where fh = m o k,^ 1 is viewed as a function on C 2 in the natural way. 

Remark. We notice that in view of (I2.12p . the inclusion map H(A,m) — > H(A) is of 
trace class. 

When defining the action of the operator P £ on H(A), we follow [23], [29], and 
perform a contour deformation in the integral representation of P e on the FBI- 
Bargmann transform side. The operator P e then receives a leading symbol given by 
Pe\a- More generally, when a E S(A$ e ,rh), with an order function m on A$ £ ~ C 2 , 
we would like to define the Weyl quantization of a = a o kt € C°°(A), acting on 
H(A). To this end, it will be convenient to introduce a globally unitary /^-Fourier 
integral operator with a complex phase, 

U = U £ :L 2 (R 2 )^H(A), U £=0 = 1, (3.54) 

depending smoothly on e, associated to the canonical transformation in f)3.47p . in 
order to reduce the quantization procedure on A to that of Weyl on T*R 2 . When 
defining the unitary operator U globally, we follow the procedure described in detail 
in [30], [25]. We can then introduce 

Op h (a) = 0(1) : H(A, fh) -> H(A), (3.55) 

defined as 

Op h (a) = Uo(aon) w (x, hD x ) o U\ (3.56) 

using the /i-Weyl quantization on R 2 . In what follows, when quantizing symbols 
defined on the IR-manifold A, we shall always use the unitary operator U to reduce 
to the standard phase space T*R 2 . 
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Remark. In the case when M is compact, the definition of the spaces H(A) and 
H(A,m) = H(A, (a^) m ) has been given in the appendix of [12], following [3TJ. We 
refer to the latter paper for a discussion of the action of /i-differential operators with 
analytic coefficients on H(A). When quantizing a symbol a £ Cg°(A), we follow [31] 
and use the Toeplitz quantization, as explained in that paper. 



The following result is a consequence of Proposition 13 . 1 1 and the arguments of Section 
5 of [16]. See also Proposition 6.3 of [T5j . 



Proposition 3.2 Assume that ( I2T29]) holds and that h K <£<h 5 ,0<5<l<K. 
Let < J< 1 be sufficiently small and let Nq > 1 be fixed. Assume that 



z e 



£ 1 f! 

C" c 



IE 



p. TP. 

3 C" 3 



C 



1,3, 



is such that dist(z, Spec(P e )) > eh N ° . We then have, in the sense of linear continuous 
operators: H(A) — > H(A,m), 



Pr 



< 



(3.57) 



4 Trace class perturbations and the vertical bo- 
undary segments 

In this section we shall be concerned with the trace of the integrals 



1 

27TZ 



— / (z-PX 1 dz, 



J 



2,4, 



->j 



and, to fix the ideas, we shall take j = 2. It will be clear that the treatment of the 
case j = 4 is similar. Furthermore, for simplicity, we shall concentrate on the case 
when M = R 2 . 

In the Hilbert space H(A), let us write 

P £ = P + 0(e) : H(A, m) -> H(A), 

where P = UP e=0 U* is selfadjoint in H(A). Here P e=0 is selfadjoint on L 2 (R 2 ). We 
shall write p £ C°°(A) for the leading symbol of P. 

Following the approach of [22], [33], we shall introduce a trace class perturbation of 
P £ in order to create a gap in the spectrum. To this end, let x £ Co°(Pt<) [0, 1]) be 
such that x(t) = 1 f° r \t\ ^ 1? an d consider the symbol 



Hp; h ) = ix 



p{p) - E 2 



-3<5 



p £ A. 



18 



We shall assume that < 8 < 1 is so small that 

for some rj > 0. Associated to k(p; h), let us introduce 

P = P + e f5 K, (4.2) 

where _?T is the Weyl quantization of k(p; h), obtained by using the h-Wejl quanti- 
zation on R 2 , as described in (13.561) . 

The hypothesis (14. ip together with the standard estimates for the operator norm 
and the trace class norm of a pseudodifferential operator on R n , given in [5|, show 
that 

||P-P|| < o(e 31 ) , ||P-P||tr < O I — I . (4.3) 



e 65 
Write next 

P £ = P £ + e 3l K. 

We shall make use of the following crude parametrix construction for z — P £ , when 
Rez = E 2 and Imz = 0(e). In doing so let us consider 

e{p,z,e) = l— y UP) :=Pe(p)+ie 3l X ( P{p) ~z E2 )- (4-4) 

z-Pe{p) V. e 35 J 

Here p £ is the leading symbol of P £ , acting on H(A). We shall first restrict the 
attention to the region where 

\ P (p)-E 2 \<0(e 31 ). (4.5) 
Here we have, with Rez = E 2 , Imz = 0(e), 

£ 38 

\Up)-A>-^y (4-6) 

Considering the usual expression for V £ e, I > 1, given by the Faa di Bruno's for- 
mula [21], we see that we have to estimate the expression 

7 -Til? li + ...+l k = l. (4-7) 

[Z ~ Pe) fJl {Z ~ Pe) 

Distinguishing the cases £j = 1 and £j > 2, and using that Vp £ = 0{1) together 
with (14. 6p . we obtain that in the region where (I4.5P holds, we have 

V e e(p,z,e)=e- 3l O e (e^). (4.8) 
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When p e A is in a bounded set away from the region where (14. 5 j) is valid, the 
estimate ( 14. 8 p improves to the following, 

Finally, when p is in a neighborhood of infinity, where A agrees with R 4 , we get 

V'e(p, £ ,;) = O f (^). (4.10) 

Introducing E(z) to be the Weyl quantization of e(p, z, e), defined using the unitary 
map (13.541) . as described in Section EJ we conclude, using (Q]) . (Ojl . (j4~9|) . (14.101) . 
that 

E(z) = O f -L) : if (A) £T(A, m), (4.11) 
while the trace class norm of the operator E(z) on H(A) satisfies 

Here p{dp) is the symplectic volume element on the IR-manifold A, so that 

p{dp) 



a 2 



2! 

where a is the complex symplectic (2,0)-form on T*M, and K C A is a sufficiently 
large fixed compact set. 

It is then clear from (Q]l . (Q|) . (jO|) . and (ODj> . that 

(z - P e )£(z) = 1 + i?, (4.13) 

where 

« = o(Jy) : 27(A) ->2f(A). (4.14) 

It follows therefore that in the region where Kez = E 2 , Imz = 0(e), the operator 
z — P £ : H(A,m) — >■ H(A) is bijective, and 

(z-P.)"^©^) :ff(A)-+fr(A,m). (4.15) 



Furthermore, 
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and writing (1 + P)" 1 = 1 - (1 + P^P, we get 

(z-P e ) 1 -E{z) = -E{z){l + R)~ 1 R. (4.16) 

It follows therefore from (I4.12p and (I4.14p that the trace class norm of the operator 
in the left hand side of (14.1 6p is 

o(^logi), (4.17) 

assuming that ( 14. ip holds. 
We conclude that 

Here we have also used that the length of 72 is 0(e). 

We shall now compare traces of the integrals of the resolvents of P £ and P e , fol- 
lowing some classical methods in non-self adjoint spectral theory. Let us also re- 
mark that such methods now also have a tradition in the theory of resonances. 
See 0, [22], [23], [32], [31], [26], [38], [42]. We shall therefore only recall the main 
features of the argument, referring to the above mentioned works for the details. 

Assume that Kez = E 2 and Imz = 0(e). An application of the resolvent identity 
shows that 

(z - Pe)^ ~(Z- Per 1 = ~{Z- PeY 1 K (z - P.)"' 

is of trace class on H(A), and using the cyclicity of the trace, we obtain, by a classical 
calculation, that 

tr Uz - PeT 1 - (z - P e ) = tr n 1 + P>)) d z K(z)^j . 
Here K(z) = e sI K (z - P £ ) • Hence, 

*{hSr- p ' rUz -^L (z -~ p ' rdz 



— [ <9 2 logdet (l + e 3S K(z- Pe)" 1 ) dz. 
2m J 72 V / 
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We shall be interested in the real part of the expression above, which is equal to 

1 

where we have set 

D(z) = det [l + K(z)\. (4.19) 



— var arg 72 D(z) 



An application of (14. 3D and (14. 15)) shows that 

h 2 j 



K(z)\\tr<0[ t -P\, (4.20) 



in the region where Re 2; = E 2 , Imz = 0(e). 

When estimating the argument variation of D(z) along 72, we shall proceed by fol- 
lowing the now well established and essentially classical complex analytic argument, 
described in detail in [22], [32], [34]. (See also [37].) In order to recall its main 
features, let us consider the holomorphic determinant D(z) in (14. 19ft in a region of 
the form 

R d = [E 2 - de 3S , E 2 + de 35 ] + i[-de 3S , de 3S ], d>0, (4.21) 
and notice that the bound f l4.20p remains valid for z G Rd- It follows that 

\D(z)\ <exp (|| K{z)\\ tT ) <exp (o{l)^j , z G R d . (4.22) 
Let now zq G Rd be such that Imzo < and |Imzo| > d\e 3S , d\ < d. We then have 



and therefore, since 
with 

it follows that 



[z -Pe)- l = o(^y.H(A)^H(A), 
Dizo)' 1 = det ((1 + X^o))- 1 ), 



[l + K{z ))- 1 = l-e^K{z -P t 



-1 

el ? 



e 3S 



Let = N(P E ,Rd, h) be the number of eigenvalues Zj, j = 1, . . . N, of P e in Rd, 
repeated according to their multiplicity. Using that (I4.20p continues to be valid in 
a region of the form Rd, with a slightly larger value of d, and combining this with 
(14.231) and Jensen's formula, we obtain that 

F 3S 

N{P e ,R d ,h) = 0(l) — . (4.24) 
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Proceeding further as in [22], [32], [34], one next considers a factorization 

JV 

D(z) = G(z) \[{z - z,), z e Rd, (4.25) 

i=i 

where G and 1/G are holomorphic in R^. An application of Cartan's lemma (or, 
alternatively, of Lemma 4.3 in [34] ) together with the maximum principle and the 
Harnack inequality allows us to show that after an arbitrarily small decrease of d > 0, 
we have 

|log|G(*)||<of^V zeR d . 

It follows then easily (see, for instance, Lemma 1.8 in [22]) that the argument vari- 
ation of G(z) along 72 is 



vararg 72 G(z) = O 



£ 36 



h 2 

Combining this with ( 14. 24ft and ( 14.251) . we obtain that 



var arg D(z) = O 



Proposition 4.1 Assume that 5 > is small enough so that hz n < e 35 , for some 
r] > 0. We have 

Re ( tr {h Jj> ~ p ^ d * -hjf- dz )) = (f l • (426) 

Combining (I4.18P and Proposition 14. 1[ we obtain that 

1 f If ( f 3 ^ hr~ 5S 

Retr— / (z-P,)' 1 dz = Retr— E(z) dz + O [— + -—log - 
2-ki J 12 2m J 12 yh 2 h 2 e 

'(4.27) 

and here the remainder in the right hand side is 

IF' 



provided that 5 > is so small that 



h < jz r. (4.2? 

loge 
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Notice that the smallness condition (14. ip is implied by (I4.28p . provided that r\ > 
in (14. ip is sufficiently small. 

Since the operator E(z) is introduced by means of the Weyl quantization on R 2 and 
the unitary map ( 13.541) . associated to the canonical transformation ( 13.4TI) . we know 
that the trace of the trace class operator E(z) : H(A) — > H(A) is given by 



Here we recall that p e has been introduced in (I4.4p . 

With (I4.27P in mind, we shall now compare the expressions 

±f [[ 1 

2niJ J2 JJ z-p e (p) 

and 



p{dp) dz (4.29) 



-Iff ^7T tidp) dz. (4.30) 

The difference between the integrals (I4.29P and (I4.30p is equal to 

1 r rr^x((p-E 2 )/e 3 A 

~ // / ~v r^^dz, (4.31) 

~lJ l2 JJ {Z-Pe){z-Pe) 



2m 

which, in view of (14. 6p . does not exceed 



which can in turn be estimated by 

0(1) fJ X {^-^ ] j{-\og\ P -E 2 \)p{dp) = 0{l)^ -\ogtdt = 0{l)e ? >hog 1 -. 

We summarize the discussion in this section in the following proposition. 

Proposition 4.2 Let E 2 < < E 4 be such that \Ej\ ~ e s , j = 2,4, where < 5 < 1 
is so small that 

I log ^ I 

When 7j is the vertical segment given by Re z = Ej, eF 3 < Imz < sF\, we have, for 
J = 2,4, 

Retr^- f (z- P £ y l dz 



2wi 



->j 



1 1 



21 



Remark. As mentioned above, the idea of using trace class perturbations to create 
a gap in the spectrum of a non-self adjoint operator has a long tradition in non- 
selfadjoint spectral theory, [23], [33]. Here we have chosen to create gaps that are 
wide enough, so that simple pseudodifferential perturbations can be employed to 
that end. The price to pay for this simplicity is that the remainder estimates that 
one obtains in the trace analysis in Proposition 14.21 are not expected to be sharp, 
and it is quite likely that finer estimates are possible to derive, at the expense of 
a greater technical investment. We hope to be able to return to this question in a 
future paper. 



5 Trace analysis near the Diophantine levels 

The purpose of this section is to understand the semiclassical behavior of the trace 
integrals 

tr ^/ (z-Pe^dz, J = 1,3, 

and when doing so, we shall concentrate on the case when j = 1. Here we recall that 
7i is given by Imz = eFi, E 2 < Rez < E±, \Ej\ ~ e 5 , j = 2,4. We shall assume 
throughout that < 5 < 1 satisfies 

h < 0{e 9 ~ s ), (5.1) 

so that Proposition 14.21 is applicable. Recall also that F% < F\. We may, and will 
assume, in the remainder of the following discussion that z G 71 satisfies 

dist(2,Spec(P £ )) > eh No , (5.2) 

for some fixed A > 1, so that by Proposition 13.21 the resolvent (z — -Pg)" 1 satisfies 
an estimate of the form (I3.57P . 



5.1 Trace integrals away from the Diophantine tori 

Following [16] , thanks to Proposition I3.1[ we shall consider a smooth partition of 
unity on the manifold A, given by 

1 = xi + X3 + A,+ + A- + + iPi,o + (5-3) 

Here < Xj £ Qj°(A-) is a cut-off function to an ^-neighborhood of Aj, j = 1,3, 
such that, in the sense of trace class operators on H(A), we have 

[P e ,Xj] = 0(h M ). (5.4) 
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supp z/v 



Rep, 



Figure 1: A schematic representation of the partition of unity (15. 3p on the IR- 
manifold A, chosen according to Proposition 13.11 Here p £ is the leading symbol of 
P e , acting on H(A). 
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The integer M = M(N, 5,5) > is fixed and can be taken arbitrarily large by 
choosing the integer N in Proposition 13.11 large enough. As observed and exploited 
in [IB] , such a choice of the cut-off Xj is possible thanks to Proposition 13.11 



The functions < ip T ,± in (I5.3P are chosen so that 



±ReP £ (p)>-^j- y (5.5) 

near the support of ift rt ±, respectively. Here we may assume that the support of 
■?/v,+ contains a neighborhood of infinity, where A agrees with T*M, and there the 
estimate (I5.5P improves to the following one, 

ReP ^ p) -W)- (5 - 6) 

We shall now describe the properties of the cut-off functions ipi± and ipio, occurring 
in ( 15. 3ft . The function < G Cq°(A) is such that near the support of ipi- we 
have 

while near supp(^j i+ ) we have 

ee 1 

lmP £ > eF 1 



0(1)' 
Finally, we have 



ee 5 . „ „ ee s 



^3 + ^ttt < ImP £ < eF 1 



0(1) - * - - - 

near the support of € Cq°(A). Continuing to follow P3], we may and will arrange 
so that in the sense of trace class operators on H(A), we have 



A[P e M = 0(h M ), - = ±,0, (5.7) 

and also, 

[P £ ,^ h ]A = 0(h M ), - = ±,0. (5.8) 

Here A is a microlocal cut-off to a region where |ReP e | < e s /0(l), and M = 
M(N, 5, 5) is an integer having the same properties as the integer in (15 .4p . 

Decomposing the operator (z — Pe) -1 according to (15. 3p . we shall first analyze 
the trace of the integral 
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We have 

(z - P,)- 1 = O (J-) : H(A) //(A, m), (5.9) 



and from [16], let us recall that 



(z - P,)- 1 W,+ = O (i=J : if (A) if (A, m), 

provided that h < 0(e 9S ). Thanks to the essentially elliptic estimate (15. 5p . it is 
possible to construct a trace class parametrix for z — P e , valid near the support of 
ip r ,+- We shall now describe briefly the main steps of this well known construction. 
Let x — Xr,+ £ C£°(A) be such that ReP e > e 5 /0(l) near the support of x an d 
X = 1 near supp ijp r ,+)- Let 

e {p,z,£) - 



Z-Pe{p)~ 

Restricting the attention to a suitable bounded region of A, we see that 

V £ e = O [e~~ & e^ , feN, 
while in a neighborhood of infinity, where A agrees with R 4 , we have 

v% = o ( - 

It follows that, on the level of operators, we have 

(z - P £ )e = X + ri, n = (J^J : if (A) H(A). 

Here and in what follows, we are quantizing the symbols on A using the Weyl quan- 
tization on R 4 , as explained in Section [3) Notice also that the trace class norm of 6q 
does not exceed 

\o (log 1 



Furthermore, with 



h 2 



-riX 

ei 



Z-Pe 

we get on the level of operators, 

h 2 



z - P £ ) (eo + e 1 ) = X + ri(l - x) + r 2 , r 2 = O ( ) : if (A) if (A). 
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and continuing in this way we obtain the symbols ej, 1 < j < L, and r^, 1 < k < 
L + 1, L G N, such that on the level of operators, 

L 

(z-P e )(e + e 1 + ... +e L )=x + ^r k (l- X ) + r L+1 . 

k=i 

Here 

r k = 0(h M *):H(A)^H(A), 

where M k — > oo as k — >■ oo. It follows, using also (15. 9ft . that modulo an expression 
whose trace class norm on H(A) can be estimated by an arbitrarily high power of h, 
provided that we take L large enough, we have 



(z - P e ) Vr,+ = (e Q + ei + ... +e L ) tp. 



r,+ - 



Estimating the trace class norm of ejxf) ri+ , j > 1 and using that the length of 71 is 
0(e s ), we obtain the following result. 

Proposition 5.1 We have 

tr ^— I (z- Pe)" 1 ij) r ± dz 
2wi y 7l 

US t4m iy±(p) Kdp) iz + log i) ■ (510) 



2ni (2iih) 

We shall next consider the trace integral 

tr-^ / (z-P £ y l ^_dz. (5.11) 
2m J 71 

Here we are no longer in the elliptic region, and in order to understand (15. lip , we 
shall proceed by means of a suitable trace class perturbation, concentrated in the 
non-elliptic region, followed by a contour deformation argument. To be precise, let 
P £ be such that P e = P e near supp(^ i _) and with ImP e < eF 3 — ee s /0(l) in the 
entire region where |ReP £ | < e s /0(l), while P e agrees with P £ outside of a slightly 
larger set of the form |ReP £ | < e s /0(l). We shall arrange, as we can, so that the 
operator P e is of the form 

P £ = P £ -ie X , (5.12) 
where \ — x(P'i h) G C£°(A) satisfies 

V e X = O e (e- e '), i>0. (5.13) 
We shall make use of the following result. 
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Lemma 5.2 Assume that 5 > is so small that h < 0(e 9S ), and recall that along 
7i, the bound dist(z, Spec(P e )) > sh N ° holds, for some fixed N > 1. We have 

1. In the region where Kez G [E^-Ei], eF\ < Ymz < 0(e s ), the operator 

z-P £ : H(A,m) ^H(A) 

is bijective, with 

z - PeY 1 = 7 : H ( A ) H(A, m). (5.14) 

/ Im z - eF 3 + ee 5 

2. It holds, in the sense of trace class operators on H(A), 

(z - P e )-Vi,- ~(z- Pe)' 1 = 0(h M ), z G 7l , (5.15) 

where M can be made arbitrarily large, by taking the integer M in (15. 7p large 
enough. 

Proof: In the region where |ReP £ | < e s /0(1), we have ImP e = 0(e) and 

1 ~ E 1 

-ImP F < F 3 



E 6 - ° 0(1) 

The statement (1) is obtained therefore by a standard application of the sharp 
Garding inequality. See also Section 5 of [16] for the details of a similar argument. 

When establishing (2), let us notice first that the expression in the left hand side 
of (I5.15P is equal to 

(Z - Pe)-\Pe ~ Pe) (z ~ P e ) _1 (5.16) 

Let ipi- be such that ipi- = 1 near supp(^ ) _) and with P £ = P £ near supp(^ i _). 
We can also arrange that (15. 7p and (I5.8P are also valid for ipi,-. Modulo 0(h M ) in 
the trace class norm, we may replace the expression in (I5.16P by 

z - p)" 1 (p £ - P E ) (l - (z - P)" 1 
whose trace class norm does not exceed 0(e~ & ) times the trace class norm of 

(l-*,-) (Z-Pe)' 1 

When estimating the latter, we may follow some arguments of [31] • When L G N, 
let 

ipi - := i> -< -01 -< . . . -< 4'l -< 4' := ipi _. 
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Here we arrange also that ( 15. 7p and (15.81) are valid for ipj, 1 < j < L. Modulo terms 
that are 0(h°°) in the trace class norm, we may write 

L 

+ {z- p e y l [p £ ^ L ](z - p^ip^l^} . . . [p e M{z - p £ )"V, 

and therefore, still modulo terms that are 0(h°°) in the trace class norm, we obtain 
that 

(l-?)(z-P e )-^ 

= (l - Pj (Z - Pj-^Miz ~ P E ) -1 [^^i-l] • • • [PeMi* ' ^)"V- 

Here we recall that for 1 < j ' < L, 

in the trace class norm, where A is a microlocal cut-off to the region where |Re P e \ < 
e 5 /0(l). Using also that 

(z - P e )-\l -A)=o(iy. H(A) -> H(A,m), 

while the trace class norm of this operator is 

W 1 



h" 1 Vel 

as well as the fact that [P £ ,ipj] = 0(h/e 25 ), we obtain the second result. □ 
An application of Lemma 15.21 shows that 

tr 2^-/ (z-PeT 1 tl)i,-dz = tv^-J ^z-P e y\i,-dz + 0(h^). (5.17) 
Here 

/ (z-PeY 1 iPi,-dz= (z-PX 1 ipi-dz, (5.18) 

where 71 is a piecewise linear contour contained in the region where Imz > eF\ 
and having the same endpoints as 71, such that the closed contour (—71) U 71 is 
the positively oriented boundary of the triangle with the third corner at the point 
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(E 2 + E±)/2 + ieF\ + ie i '. When understanding the integral in the right hand side 
of (I5.18p . it will be convenient to decompose the function further, according to 
the values taken by ReP e , so that we write 

= -0i-,2 + ^i-,4 + ^i,-- (5.19) 

Here the functions "0i,-,j are supported in the regions where |ReP e - E s \ < e 3S /0(l), 
j = 2,4, respectively, while the support of ipi- stays away from these regions. 

We shall consider first the integral 

z-PX 1 A,-dz. (5.20) 



71 



In the support of ipi-, we have, 

z-Ps(p) >^y, 2671, (5-21) 

which is an essentially elliptic estimate. Using ( I5.12p . (I5.13p . and (I5.2ip we then 
obtain that 

[(z - P e (p)r%,-) = e^O (e-™) , t G N. 
Also, the trace class norm of the corresponding Weyl quantization does not exceed 

Here we may replace the function by another cut-off with a slightly larger sup- 
port. It is therefore clear that the integral in (I5.20p can be understood by construct- 
ing an /z-pseudodifferential parametrix for z — P £ , valid near the support of ipi-, by 
following the same method as in the proof of Proposition 14.11 We therefore conclude 
that the trace integral 

trJLJ^ (z-P e ) '^-dz 

is equal to 



^(p)p(dp)dz + o[^—\. (5.22) 



2vri(27r/i) 2 J % J J z-p £ {pY' ™ I h 2 

Here p £ is the leading symbol of P e , where we know that p £ —p £ = 0{e) is supported in 
a region where |ReP e | < e s /0{l), away from the set where ImP e < eF 3 — ee s fO(l). 

According to ( 15. 18ft and (15.191) . it remains to consider the integrals 

[(z-P^^jdz, j = 2,4. (5.23) 
J ~n 
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Here an application of the bound (I5.14p in Lemma 15.21 together with the fact that 
the trace class norm of on H(A) for j = 2,4 is 



O 



£ 38 



h 2 I ' 

shows that the trace class norm of the expressions in (I5.23|) does not exceed 

( e 3l \ f 1 1 e z ~ 5 1 

Combining flBTTTD . floTTSD . (ICTjj) . f l5T22l) . and ( K2M we conclude that the trace 

tr-L f (z-P,)- 1 $ it _dz (5.25) 
27" J 



is given by 



1 l /■/■/* 1 ~ e 35 1 

^,_(p) //(dp) + 0(1)— log -, (5.26) 



27TZ (27r/l) 2 J^J J Z — p e (p) ' h 2 E 

provided that the lower bound (15.11) is strengthened to the following one, 

h<e x ™\og-. (5.27) 

Here the integral 

4>i,-(p) Kdp) dz 



71 



Z ~ Pe(p) 



can be replaced by the integral 

1 



71 



4>i-(p) p{dp) dz, 



Z-Ps(p) 

at the expense of an additional error not exceeding 



The result obtained so far is typical of the behavior of the trace integrals in question 
in the non-elliptic region away from the tori Aj, and so we state it in the following 
proposition. 
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Proposition 5.3 Assume that < 8 < 1 is such that h < e 105 \og^, and recall that 

^P.<eF,-^- y 
near supp (V'i,-); where ipi- satisfies (15 .7p , (15. 8p . VFe /jai>e 

US " (<w * + ° (£) log (528) 



2m' (2mh) 



It is now easy to extend the result of Proposition 15.31 to the trace integral 



tr ■ 

2m 



hi J t*-^ 1 dz - 



Indeed, the argument leading to the asymptotic result (1 5 . 2 8 1) remains valid also when 
considering the trace integrals 



7i 



and 



2vn J 



'a 

where we may also recall that F% < F±. 
When considering the expression 

/ (z - PeY 1 ipi,+ dz, 

J 71 

we introduce a new trace class perturbation P £ of P e , such that P e = P £ near 
supp(^j i+ ), and with 

ee 1 

ImP e > eFi 



0(1) 



in the entire region where |ReP e | < ^/^(l), and such that P e agrees with P e further 
away from this set. The natural analogue of Lemma 15.21 continues to be valid for 
(z — P £ ) _1 , and when studying the trace of 



/ (z-p £ y l ^ h+ dz, 

J 7! 
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we can therefore deform the contour 7 X downwards and introduce the decomposition 
similar to (I5.19p . so that the estimate f !5.2ip holds for z — P £ along the deformed 
contour, near the support of ^ )+ . 

The following is the main result of this subsection. 

Proposition 5.4 Assume that E 2 < < E±, \Ej\ ~ e 5 , j = 2,4, where < 5 < 1 is 
so small that h < e 10S \og(l / e) . Let jj, j = 1,3 be the horizontal segment given by 
E 2 <Kez < E4, Imz = eFj. Let finally < Xj £ Co°(A) be a cut-off function to an 
e 5 -neighborhood of Aj, j = 1,3, enjoying the commutator property (15.41) . We have 

tr-L / {z-PeY 1 (l-Xj)dz 

= h (2^F i II ^U) {1 " Xj(p)) Kdp) dz + ° (^) log ? 3 = h 3 ' 

5.2 The Birkhoff normal form and trace integrals near the 
tori 

In this subsection, we shall complete the trace analysis near the Diophantine levels 
by studying the integrals 




say, when j = 1 . 

When 2 G 71, let us consider the equation 

[z — P £ ) u = v, 116 H(A, m), 

so that 

- Pe) XlU = XlV + [Xl, P e ]u. 

Applying the operator U := U\ of Proposition 13.11 we get 

[z - P[ N) - R N +i,i) Uxiu = Ux\v + T N u, 

with the trace class norm of T N on H(A) being 0(h M ), where M can be taken 
arbitrarily large, by taking the integer N in Proposition [3J] large enough. Using the 
fact that i?jv+i,i(a;, £, e; h) = O ((h,e,£) N+1 ) and modifying T/v slightly, we get 
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It is therefore clear, in view of the fact that e > h , K > 1 fixed and since the bound 
(B holds, that 



Xi (z - PeY 1 = v[z- P[ N) ) ~ U X i + T liN , (5.29) 

where T^jv has the same trace class norm bound as T/y. Here 1/ is a microlocal 
inverse of U. Therefore, 

tr — / (z- P,)" 1 

= tr— / fz-P^V 1 f/ Xl F^ + C»(/i M ). (5.30) 

Here we may recall, following [IB] , that UxiV = x{hD x / 'e s ), where x £ Co°(R 2 ) 
is a standard cut-off to a neighborhood of £ = 0. Let us write in what follows 
Xe(C) = X^/e 1 )- 

Now the eigenvalues of the translation invariant operator P[ N \hD x ,e] h), acting on 
Lg(T 2 ) are given by 



fx(k) := Pr> ( h I k - -f I - ^ I , fc G Z', (5.31) 
and by (15.301) . we conclude that 
tr-^ / {z-P £ )- l Xl dz 



We have 



where 



Ai(*) =p[h[k-^\-^-,e)+ O(h), (5.33) 



p(^)=p(0 + ^>(0 + <^ 2 ) 

is the leading symbol of P[ N \hD x ,e; h). We would like to compare the expression 
in the right hand side of (I5.32p with the integral 



1 1 



J jj \z-p{^e)Y 1 Xe{Odxdidz, (5.34) 



2ni (2vr/i) 2 

the integration in the (x,£) variables being carried out over T*T 2 . 
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Integrating out the x- variables in (I5.34p . we shall first compare the expressions 

[ [ (z-p(^e)y 1 Xe(Z)d£dz 



and 



£'i(-K fc (»-T)-i-)r-M*-s-s)*- 



fcez^ •''ft 
When £ e R 2 is such that 



2 



for some fcGZ , let us write 
Let us consider 

/ f {(z-p^eT'-iz-pd^eT^XeiO^dz. (5.35) 

Let a, b be the endpoints of 71. Then with suitable branches of the logarithm, we 
have 



(log(6-p(£,e))-log(&-p([£],e))) 

- (log (a - £ )) - log (a - p([£], e))) . (5.36) 



In general, for z, w e C, we have 

[ w 1 

\ogz-\ogw = J -d(, 

where the choice of curve joining w and z depends on the choices of branches of log z, 
log w. If we have the same branch then 



Co 


z ■ 




10 




min( \z 


> 




1) 



logz — logu>| < — — : — pr. (5.37) 



If the branch cut passes between z and w, we have to add a constant. In the case of 
f)5.36p . this happens precisely when p[£,e) and p( [£],£:) are on the opposite sides of 
71. Now let us concentrate on one of the terms in (I5.36p . say 

log (a - p(£, e)) - log (a - p([£], e)) . (5.38) 
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If 

\Rea-Rep(£,e)\<C h, (5.39) 

for a suitable fixed constant Co > 0, we estimate the two terms separately and get 
that the contribution to ( 15.351) in this case is 

/ x B (Olog(a-p(£,e))d£<0(l) / * -logtdt = 0(l)h log \. 
Je(c q ) Jo n 

Here E{Cq) C R 2 is the set of all £ G R 2 such that (15.391) holds. If we assume that 
a has been chosen so that for all £ G R 2 , 

|a-?([&e)l>on), (5-40) 

then we get the same estimate for 

Xe (e) log (a -£([£], £ )) dC 

In the region where |Rea — Rep(£, e)| > Co/i, let us first assume that we have the 
same branches of log(a — p(£, e)) and log(a — e)). Then by (15.371) . 

log( a -^,e))-log(a -?([{],,)) = (i) 5 _^_ J , 

and the corresponding contribution to (I5.35P is 

1 1 . _ (. . 1 



Jh t V " 

It remains to estimate the integral of the extra contributions ±2tri, to log(a— e)) — 
log(a— e)) from points£for which p(£,e) and p( [£],£) are on the opposite sides of 
7i. But the Lebesgue measure of the set of such points is 0(h), so the corresponding 
contribution to the integral is 0(h). 

Summing up our estimates and computations, we see that the expression (I5.35P is 
0(/ilogr). Arguing similarly and using (I5.33p . we obtain that 

/ / ({z-miz)T X -{z-KW X )xm)didz = O{h\0g\). (5.41) 

J-R? J ix n 

Finally, we find that also, 

/ / (z-p([^e)r 1 (Xe(0-Xe(md^dz = O(h\ogh. (5.42) 

We summarize the result of this subsection in the following proposition. Here we also 
use that the integral over T*T 2 in (15.341) can be transformed into the corresponding 
integral over A by means of the canonical transformation associated to the operator 
U. 
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Proposition 5.5 Assume that E 2 < < E A , \Ej\ ~ e 5 , j = 2,4, where < 5 < 1 is 
so small that h < e 105 log(l / e) . Let jj, j — 1,3 6e i/je horizontal segment given by 
E 2 < Kez < E±, Imz = eFj. Let finally < Xj G Co°(^-) ^ e a cut-off function to an 
e s -neighborhood ofAj, j = 1,3, enjoying the commutator property (15.41) . We /iave 

= in(2^UI J^pKp)^ ^ dp) dZ + ° (1) bg ? J = lj 3 - 
5.3 End of the proof 



Combining Propositions I4.2[ I5.4[ and I5.5[ we obtain that the number of eigenvalues 
of P e in the rectangle R in (I2.33P is equal to 

Here fl(e, R) = C A and p e G C°°(A) is the leading symbol of P e , acting on 

H(A). 

Let us now recall the C°° canonical transformation 

k : neigh(p _1 (0), T*M) -> neigh(p _1 (0), A), 

introduced in (13.471) . so that k is (9(e)-close to the identity in the C°°-sense, for 
each fixed T > T , T > large enough. An application of (13.481) shows that 

pMp)) = P(P) +ie(q- H P G) (p) + 0(e 2 ), p G T*M. (5.44) 
In the compact case, we obtain the same expression for the transformed symbol. 

It follows from (I5.44p and the properties of the function G, recalled in Section [3l that 
the set Q(e, R) C A is C(e)-close to the set 

n([E 2 ,E 4 ]):= [j Q(E) c T*M, 

E 2 <E<Ea 

introduced in (12.351) . Therefore, 

fi(dp)= // dxd£ + 0(e). (5.45) 

Here T > is large enough fixed. Combining (I5.43P and (I5.45p . we complete the 
proof of Theorem 12. 1[ 
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6 Weyl asymptotics in the periodic case 



In this section, we shall explain how the results and methods of the work 
combined with the methods of the present work, can be used to obtain an analog 
of Theorem 12.11 in the case when instead of the complete integrability assumptions, 
we assume that Hamilton flow of p is periodic. It turns out that the analysis of the 
periodic case will proceed in full analogy with the previously analyzed completely 
integrable case. The following discussion will therefore be somewhat brief. 

In order to fix the ideas, throughout this section, we shall consider the case when 
M = R 2 . LetP £ be an operator satisfying all the assumptions made in subsection 2.1, 
and in particular, (j2.12jl . Assume that for E e neigh(0,R), the following condition 
holds, 

The ifp-flow is periodic on p~ 1 (E) D R 4 with (6.1) 
period T(E) > depending analytically on E. 

As in (EOT]) , we set 

I rT(E) 

(q) = j^)J qoexp(tH p )dt on p~ 1 (E) H R 4 , (6.2) 

and notice that the functions p and (q) are in involution, so that H p (q) = 0. Similarly 
to ( pIUD , it is established in [12] that, 



-Im (Spec(P e ) n {z; |Rez| < 5}) C 



min (q) — o(i), max (q) + o(l) 
P -!(o)nR 4 p-!(0)nR 4 



(6.3) 



as £, h, 5 —7- 0. 
Let 

Fj e [ min (g), max (q)), j = 1,3, F 3 < F u 
V!(o)nR 4 P -!(o)nR 4 

and let us introduce the associated level sets 

A , F . = {pe R 4 ; p{p) = 0, (q){p) = Fj}, j = 1,3. (6.4) 
As in [12], we shall work under the general assumption that for j = 1,3, 

T(0) is the minimal period of every ifp-trajectory in A 0t p.. (6-5) 
We shall furthermore assume that 

dp, d{q) are linearly independent at every point of A ^ , j = 1, 3. (6.6) 
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It follows that each connected component of the level set A 0> f is a two-dimensional 
Lagrangian torus. For simplicity, we shall assume that the sets A 0i f are both con- 
nected, j = 1,3. We can then make a real analytic canonical transformation 

Kj : neigh(A , F:P R 4 ) -»• neigh(£ = 0,T*T 2 ), 

given by the action-angle coordinates near Aq^ , so that when expressed in terms of 
the coordinates x and £ on the torus side, we have ponj 1 = Pj(^i), (q)oKj 1 = (<£,)(£)• 

Following the analysis carried out in Section 4 of [T2], we shall now state the following 
result, which can be viewed as an analog of Proposition 13.11 in the present periodic 
situation. 

Proposition 6.1 Let us keep all the general assumptions of subsection 2.1 and make 
furthermore the assumptions (16. ip . ( I6.5p . and ( 16. 6p . Assume that e = 0{h 5 ), < 
8 < 1, is such h/e < 5q, for some < 5q sufficiently small. There exists a smooth 
IR-manifold A C C 4 and smooth Lagrangian tori Ai and A3 C A, such that when 
p G A is away from a small but fixed neighborhood of Ai U A 3 in A , we have 

\ReP e (p)\>-^-r (6.7) 

or 

\lraP £ {p)-eF x \>-^- and \lmP e (p) - eF 3 \ > (6.8) 

The manifold A is an 0{e + ^-perturbation ofH 4 in the C 00 -sense, and it agrees 
with R 4 outside of a neighborhood o/p _1 (0) PI R 4 . We have 

P £ = 0(1) : H(A,m) ->■ H(A). 

When j = 1,3, there exists an elliptic h-Fourier integral operator 

U j ^O(l):H(A)^L 2 e (T 2 ), 

such that microlocally near Aj, we have 

W = 

Here Pj = Pj{hD x ,e, J; h) is an operator acting on Lg(T 2 ) with the symbol 
P 3 (c, e, - £ - h^j ~ Pj (d) + e h k r hk (t, e, ^ , |£| < -± 



0(1) 



where 



s 



and 

r 
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From [12], we also infer that the natural analog of Proposition 13.21 is valid, when the 
spectral parameter z belongs to the rectangle 



"11" 






+ is 


~C' c_ 



3 c 3 c 



j 



1,3, 



for a sufficiently large constant C > 0, and is such that dist(z, Spec(P e )) > eh N °, for 
some Nq > 1 fixed. 



Let 



R={E 2 ,E 4 )+i£{F 3 ,F l ), 



(6.9) 



where £2 < < -E4 are independent of h, with \Ej\ < 1/C, j = 2,4, for C > 
sufficiently large. We decompose dR as in (12.411) . and notice that the analysis of 
Section @] applies to the traces of the integrals 



1 

2ni 



— / (z-PX l dz, 



J 



2,4, 



as it stands, so that the statement of Proposition 14.21 continues to hold in the present 
situation. For future reference, we state it as the following result. In its formulation, 
the notation fi(dp) stands for the symplectic volume element on A, and p £ is the 
leading symbol of P £ , viewed as an unbounded operator on H(A). 

Proposition 6.2 Let E 2 < < £4 be independent of h and such that \Ej\ < 1/C, 
for C > sufficiently large, j = 1,2. Assume that e = 0{h 5 ), < 5 < 1, is such 
that h/e < 5q, for some 5q > small enough. Assume also that < 5 < 1 is such 
that h < e 95 . When ■jj is the vertical segment given by Re z = Ej, eF 3 < Imz < eF 1 , 
we have, for j = 2, 4, 



Retr 



1 

2~iu 



P. 



Re 



_1 dz 
1 



2m (2nh) 2 



->j 



z _ l p ( , fJ.(dp) dz + (e 3 hog X - 



When analyzing the trace integral 



tr-^ / {z-P £ y l dz, j = l,3, 
2m 



say, with j = 1, we use Proposition 16.11 to introduce the following smooth partition 
of unity on the manifold A, similar to (15.31) . 

1 = Xl + X3 + 1p r ,+ + A - + fa - + V'i.O + - (6.10) 
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Here < Xj C^°(A) is supported in a small flow- invariant neighborhood of A,,-, 
where P e is intertwined with Pj, according to Proposition 16. 1[ and x = 1 near Aj, 
j = 1,3. As in [13], we choose Xj so that in the sense of trace class operators on 
if (A), we have 

[P„Xj] = 0(h°°). (6.11) 
The functions < Vv,± £ Cb°(A) are such that 

±ReP £ (p)>^, (6.12) 

near the support of Vv,±> respectively. Here, as in Section El we may assume that 
the support of -?/v i+ is unbounded, and near infinity, the bound (I6.12p improves to 

RePe{p) ~W)- (6 - 13) 

We now come to describe the properties of the functions and z/^o in ( I6.10p . 
These non- negative functions in C£°(A) are supported in regions invariant under the 
ifp-flow, and the estimate 

lmP £ <eF 3 -^ 
holds near suppipi-. Similarly, near supp^j i+ we have 

lmP £ >eF 1 + ^ 1 

and finally, the bound 

£j p 3 + £.<i mPe < £Fl _£_ 

is valid in a neighborhood of the support of ^«,o- As in Section |5l we may and will 
arrange so that in the sense of trace class operators on H(A), we have 

A{P £ ,^.} = 0(h°°), [P £ ,i> i ,,]A = 0(h°°), - = ±,0. (6.14) 

Here A is a microlocal cut-off to a region where |ReP £ | < 1/(9(1). 
The analysis of the trace integrals 

tr^r / (z-P £ y l iP rj± dz (6.15) 



7i 



proceeds exactly as in the proof of Proposition I5.1[ thanks to the elliptic estimates 
(I6.12p . (I6.13p . and as there, we find that 



tr 

2 



— / (Z - P £ ) 1 1p ri± dz 

:JJl7^kur) AAp)Mdp)dz + {ji)- (616) 



2m (2nhY 
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When understanding the trace 



- f (z-Pj-ty-dz, 



2m 

we continue to follow the analysis of Section and introduce an auxiliary trace class 
perturbation of P E , concentrated in a region of the form |ReP £ | < 1/0(1), similar to 
(j5.12p . The arguments of Section \5\ apply then as the stand, and we get the following 
direct analog of Proposition 15.41 

Proposition 6.3 Assume that E 2 < < E 4 , \Ej\ < 1/0(1), j = 2,4, and let 
< 5 < 1 be so small that h < e 12S . Let jj, j — 1, 3 be the horizontal segment given 
by E2 < Re 2 < E4, Imz = eFj. Let finally < Xj Co°(A-) &e a cut-off function 
to an e 5 -neighborhood of Aj, j = 1,3, enjoying the commutator property (16.1 ip . We 
/law e 

- him* I IS ^tm (1 - *<"» * + (7?) 108 ? i =1 > 3 - 

Here we continue to assume that h<^e< 0(h 5 ), 5>0. 

An inspection, using the normal forms near the tori Aj, j = 1,3, described in Propo- 
sition 16.11 shows next that the result of Proposition 15.51 remains valid in the present 
situation. Combining this observation with Propositions 16.21 and 16.31 we conclude 
that the number of eigenvalues of P £ in the rectangle R in (16. 9p is given by 

l R (p £ (p)) n(dp) + ±=0 (e 3 hog- 



(2tt/i) 2 J J n ^ F '^" ^ h2 

Recalling the construction of the IR-manifold A, described in detail in [12], we may 
summarize the discussion in this section in the following result, analogous to Theorem 
12.11 in the periodic case. 

Theorem 6.4 Let P £ satisfy the general assumptions of subsection 2.1, in particular 
( 12.12p . and make further the assumption ( 16. ip . Let 

Fa e [ min iq), max (q)], j = 1,3, F 3 < F x , 

be such that the assumptions (16.51) and (16. 6p are satisfied. Furthermore, assume 
that the level sets Ao,Fj in (16.41) are connected, j = 1,3. Assume that e = 0{h 5 ), 
< 8 < 1, is such that h/e <C 1. Let C > be sufficiently large and let E2 < < E4 
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be independent of h with \Ej\ < 1/C, j = 2,4. Assume finally that 5 G (0, 1) is so 
small that h < e 12S . Then the number of eigenvalues of P e in the rectangle 

R={E 2 ,E 4 ) + ie{F 3 ,F 1 ), 

counted with the algebraic multiplicities, is equal to 

7 Complex perturbations and the damped wave 
equation on a surface of revolution 

The purpose of this final section is to illustrate how Theorems 12.11 and 12.21 apply in 
the case when M is an analytic surface of revolution in R 3 , and 

P £ = -h 2 A + ieq, (7.1) 

where A is the Laplace-Beltrami operator and q is an analytic function on M. When 
doing so, we shall restrict the attention to the same class of surfaces of revolution as 
in [16], [Ej, and begin by recalling the assumptions made on M in these previous 
works. 

Let us normalize M so that the X3-axis is its axis of revolution, and parametrize it 
by the cylinder [0, L] x S\ L > 0, 

[0,L] x S 1 3 M) ^ (f(s) cos9, f(s) sine, g(s)). (7.2) 

Here the parameter s G [0, L] is the arclength along the meridians, so that {f'(s)) 2 + 
(g'(s)) 2 = 1. The functions / and g are assumed to be real analytic on [0, L], and 
we shall assume that for each k G N, 

/( 2fc )(0) = f (2k \L) = 0, 

and that /'(0) = 1, f'(L) = — 1. As we recalled in [IB], these assumptions guarantee 
the regularity of M at the poles. We assume furthermore that M is a simple surface 
of revolution, in the sense that f(s) > on (0, L) has precisely one critical point 
so G (0, L), which is a non-degenerate maximum, so that /"(so) < 0. To fix the ideas, 
assume that f(s ) = 1. Associated to s we have the equatorial geodesic C M, given 
by s = so, 9 G S 1 . 

Writing 

T* (M\{(0, 0, g(0)), (0, 0, g(L))}) ~ T* ((0, L) x S 1 ) , 
and using (]7.2p . we see that the leading symbol of Pq = —h 2 A on M is given by 

p( s ,e,a,e*)=o- 2 + ^- (7.3) 
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Here a and 9* are the dual variables to s and 6, respectively. Since the function p in 
(17. 3p does not depend on 9, it follows that p and 9* are in involution, and we recover 
the well-known fact that the geodesic flow on M is completely integrable. 

Let E > and |F| < E 1 ^ 2 , F ^ 0. Then the set 

A EiF :p = E, 9*=F, 

is an analytic Lagrangian torus contained inside the real energy surface p~ 1 (E). Ge- 
ometrically, the torus Ae,f consists of geodesies contained between and intersecting 
tangentially the parallels s±(E,F) on M defined by the equation 

f(s±(E,F)) 



E 1 / 2 ' 



For F — 0, the parallels reduce to the two poles and we obtain a torus consisting 
of a family of meridians. The case \F\ = E 1 / 2 is degenerate and corresponds to the 
equator s = Sq, traversed with the two different orientations. Writing A a := A l a , we 
get a decomposition as in (I2.16p . 

p- 1 (l)f]T*M= |J A , 

with J= [-1,1], S = {±1}. 

In [16], we have derived an explicit expression for the rotation number u(A a ) of the 
torus A a C J9 _1 (l), 7^ a G (—1, 1), given by 

a r +(a) 1 / a 2 \- l/2 



oo(A a ) = - l_ ds, f(s ± (a)) = \a\. (7.4) 

7T y,_ (0 ) r(s) v r{s)J 

We are going to assume that the analytic function (—1,1) 3 a H- cu(A a ) is not 
identically constant. 

Let q = q(s, 9) be a real-valued analytic function on M, which we shall view as a 
function on T*M. Associated to each a G J, we introduce the compact interval 
Qoo{A a ) C R defined as in (12.221) . We also define an analytic function 

(-l,l)9a^(g)(A a ), 

obtained by averaging q over the invariant tori A a . Let us assume that the function 
a i — y (q)(A a ) is not identically constant. 

Example. Assume that q = q(s) depends on s only. Then it was shown in [16] that 
for all 7^ a G (—1, 1), we have 

Q OQ (A a ) = {(q}(A a )} 1 



16 



where 

J(q a) /-Ma) ff s \ 

mK) = K^ry J{M = L ia) ^ iPisTW^ ^ /(s±(a)) = |a| - 

Generalizing the example above, an analytic family of the form 

q v (s,9) = q (s) + 7731(5,0), q = q, < 77 < 1, 
was also considered in [16], and it was shown that the set 

U Qoo,„(A«), s = {±1}, 

Aa6w-!(Q)U5 

has a small measure compared with that of {{q v )(A a ); a G (—1, 1)}, and that conse- 
quently, there exists a rich set of values F G U ae jQ 00ir; (A a ) satisfying the assumptions 
(123711 . (EJSiD, Pi , and (E3SD for g^. Here Q^iK) is the range of the limit of 
(Qti)t, as T — )• oo, along A a . 

The following result is a consequence of Theorem 12.11 

Proposition 7.1 Assume that M C R 3 zs a simple analytic surface of revolution 
given by the parametrization (\7.2\i , for which the rotation number a; (A), A = A a , 
defined in (I7.4p . is not identically constant. Consider an operator of the form P e = 
—h 2 A + ieq, with q = q v as above, such that the analytic function (—1, 1) 9 s 4 
(q)(A), given by the torus averages of q, is not identically constant and extends 
continuously to [—1,1]. There exists a subset 8 V C Uagj^ 00 ^) of measure fi(rj) 
tending to zero with rj, such that the conclusion of Theorem \2.1\ holds uniformly for 
£ Uagj^ 00 ^) \ ^ r )' an d 9^ ves the number of eigenvalues of P e in the region 

[E 2 , E 4 ] + is [F 3 , Fx] , E 2 < 1< E 4 , \E S - 1\ ~ el 

Remark. A more precise description of the complement of the set S v in Proposition 
7.1, given in terms of (a, j3, 7)-good values, can be found in Section 7 of [T6] . 

Remark. Let us remark finally that when discussing the operator P £ given by ( 17. lft , 
it would have been possible to allow the analytic function q on M to depend holo- 
morphically on the spectral parameter z G neigh(l, C), with q real-valued for z real 
- see also the discussion in Section 6 of [13] for a similar observation in the periodic 
case. Such an extension is motivated by the problem of studying spectral asymp- 
totics for the damped wave equation with an analytic damping coefficient, considered 
on the analytic surface of revolution M. 

We shall now apply this remark to the situation of the damped wave equation 
(-A + 2a(x)d t + df)v{t, x) = 0, on R x M, 
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where a(x) is analytic and real-valued. An important role is played here by the 
stationary solutions e lTt u(x), u ^ and the corresponding eigenfrequencies r G C, 
given by the equation 

(-A + 2ia(x)r - t 2 )u(x) = 0, (7.5) 

and we are interested in the asymptotic distribution of the eigenfrequencies r. We 
know that the large eigenfrequencies are confined to a band along the real axis (see 
[20] . [33] for such results under more general assumptions), and that the real parts 
obey the same Weyl law as for the corresponding Helmholtz equation (—A — r 2 )u = 0. 
Less is known about the asymptotic distribution of the imaginary parts, and here 
we can apply Proposition 17.11 and the subsequent remark. 

The set of eigenfrequencies is symmetric under reflection in the imaginary axis, so 
we can concentrate on the case when Rer > 1. We make a semiclassical reduction 
by putting r = w/h, 0</i<l, Rew~l and get 

{-h 2 A + 2ihwa(x) -w 2 )u = 0, (7.6) 

and we can apply our results with z = w 2 , e = h, q = 2\fza, and with P = —h 2 A as 
the unperturbed operator with leading symbol p(x,£) = £ 2 . It may now also be useful 
to recall the following general bounds on the imaginary part of an eigenfrequency w 
in (USD, 



lim inf (aW — o(l) < Imw < h \ lim sup (a)r + o(l) 
r-KJop-i(i) J yr^oo p _i (1) 

This result was obtained in [20J — see also (I2.24p for the present completely integrable 
case. 

The function (a)^ = lim T ^ 00 (a)r is homogeneous of degree in £, thanks to the 
homogeneity properties of the H p -&ow. This means that the set E 2 < Kew < E4, 
hF 3 < Imw < hF x corresponds to the set G T*M; E 2 < £ 2 < E%, F 3 < 

(a) 00 < Fi} and we notice that the conditions imposed on Fj (i.e. on the properties 
of the corresponding torus, where p = E, (a)^ = Fj) in Proposition 17.11 in the case 
when q = a, are independent of the real energy p = E. Applying Proposition 17.11 
and the subsequent remark, we get 

Theorem 7.2 Consider the stationary damped wave equation in the equivalent forms 
( 17. 5ft and (I7.6| . Assume that the assumptions of Proposition I7TT1 are fulfilled with 
q = a. Then uniformly for Fi,F 3 G UagJ < 3°°(^') \ ^v> ^ e num ber of eigenfrequen- 
cies w of (17.61) in the region 

[E 2 , £ 4 ] + ih[F 3 , F 1 ], E 2 <1< E i} \Ej -l\~h*, 
is equal to (2nh)~ 2 times 



vol{(s,0 G T*M; E 2 <^< E 2 , Q oc (A ( x, ±- \) C [Fg,^]} + 0(h sd \n-). 
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More generally, the number of eigenfrequencies w of (17. 6p in the region 
{E 2 ,E A } + z/i[F 3 ,Fi], E 2 < 1< £ 4 , 1^1 ~ 1, 

is egim/ to 

^vol {(x,£) G T*M; £ 2 < £ 2 < £ 2 , Q oc (A(x,4)) C [F 3 , F 1 ]}) + 0(^" 2 ) 



(2tt/i) 2 v 2 - - 4 ' v v ier 

TTie number of eigenfrequencies r of (17.5)) in £/ie region [E 2 ,E A ] + z'[F 3 ,Fi], where 
E 2 < £?4, Ej ~ l/h^> 1, is equal to 

n2 ^ f 2 ^ 7-i2 -o f\ f„ ^ \ \ r- \ T? T? 11 1 i /r*/i,25-2> 



vol {(x, G T*M; £ 2 < £ 2 < ^i, Qoo(A x, ) C [F 3 , Fx]} + 



(2tt) 2 V > - > - "4, ^» |^| 
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